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Abstract 

Considering a singularly perturbed convection-diffusion problem, we present an anal- 
ysis for a superconvergence result using pointwise interpolation of GauB-Lobatto type for 
higher-order streamline diffusion FEM. We show a useful connection between two differ- 
ent types of interpolation, namely a vertex-edge-cell interpolant and a pointwise interpolant. 
Moreover, different postprocessing operators are analysed and applied to model problems. 
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1 Introduction 

Consider the convection-diffusion problem given by 

-£Au-bu x + cu = /, in£2=(0, l) 2 (1.1a) 

m = 0, ondQ. (Lib) 

where b > /3 > 0, c + \b x > y > and < £ <C 1. Note that condition on c can always be 

fulfilled by a transformation v = exp(TOc)w for a suitably chosen K. 

In Q it was shown that for bilinear elements and a standard Galerkin method its solution u N 

fulfils on a piecewise uniform Shishkin mesh with N mesh cells in each coordinate direction the 

estimates 



1 1 1« - m" 1 1 | e < CAT 1 IniV, and | \u l -i/ y \ \ \ £ < C(N^ l \nNY, 

where u 1 is the standard piecewise bilinear interpolant of u and the energy-norm 1 1 1 ■ 1 1 1 £ is defined 
as 

IIMII £ = (£||V^ + 7Nlo) 1/2 - 
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Here and throughout the paper we denote by || ■ ||o the standard L2-norm on Q. and by C a generic 
constant independent of e and N. 

A property like the one above is called supercloseness and by a simple postprocessing routine 
P JH it can be exploited as interpolantwise superconvergence (for the naming convention see 
OH) in 

\\u-Pu N \\\ e <C{N^ l \nN) 2 . 

In J6l a similar result was obtained for a streamline diffusion method ifTTTl under some restric- 
tions on the stabilisation parameters. For higher order methods using (^-elements with p > 1 
so far only for the streamline diffusion method supercloseness results are known. In [3| it was 
proven that 

W\rfu-u N \\\ £ <C(N- l \nN) p+l / 2 \nN 

holds for the streamline diffusion solution u N in the case of <2 p -elements on a suitable piecewise 
uniform Shishkin mesh and conditions on the stabilisation parameters. The interpolant itp is a 
so called vertex-edge-cell interpolant ll9l[T5l. 

In El the higher order case was investigated numerically and three unproven phenomena were 
shown. First, there seems to be a supercloseness result for pointwise interpolation w.r.t. GauB- 
Lobatto points. Second, the supercloseness order is actually p + 1 and not only p + 1/2. And 
finally, these results do also hold for standard, unstabilised Galerkin FEM. In the present paper 
we prove the first of these numerical results. 

The theoretical results presented in this paper require regularity of the exact solution u and 
use solution decompositions given e.g. in H121ll3[[T6l . We assume the necessary compatibility 
conditions and smoothness of the data to be fulfilled. 

The paper is organised as follows. In Section |2] we introduce a class of layer-adapted meshes 
used in discretising the differential equation. Moreover, a solution decomposition exploited later 
in the analysis is presented. Section |3] contains the definition and analysis of some properties 
of two different interpolation operators. In Section |4] we prove supercloseness properties of our 
two interpolation operators and in the following Section [5] superconvergent numerical solutions 
are generated by postprocessing. Finally, Section |6]contains a numerical example verifying the 
theoretical results. 



2 The Mesh and a Solution Decomposition 

We define the underlying mesh as a member of the general class of S-type meshes [18]. Let 
the number N > 4 of mesh cells in each direction be divisible by 4 and a user-chosen positive 
parameter o > be given. Assume 

In practice, this assumption is no restriction as otherwise N would be exponentially large com- 
pared with e. In the latter case the analysis could be done in a standard way. We now define 
mesh transition parameters by 

X x :=^-\nN<- and Xy,:= o^/eXnN < -. (2.2) 
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Figure 1 : Shishkin mesh T 8 of £2, the bold lines indicate the boundaries of the subdomains. 
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The domain Q. is dissected by a tensor product mesh according to 

i = 0,...,N/2, 

i = N/2,...,N, 

NJ , j = 0,...,N/4, 

(l-2A y )(^-l) + i, j=N/4,...,3N/4, 
l-ay/e^(2-^j, j = 3N/4,...,N, 



Xi :- 



yj'- z 



1-2(1-A X )(1 
'ay/iff^- 



and the final mesh r w is constructed by drawing lines parallel to the coordinate axes through 
these mesh points. The function <j> is a monotonically increasing, mesh-generating function 
satisfying 0(O)=O and 0(l/2)=lniV. Given an arbitrary function fulfilling these conditions, 
an S-type mesh is defined and the domain Q. is divided into the subdomains Qn, ^12, ^21 , and 
^22 as shown in Figured] Related to the mesh-generating function 0, we define the mesh- 
characterising function y = e~^ whose derivative yields information on the approximation 
quality of the mesh, usually expressed in terms of max | y/\ where the maximum is taken over 
t G [0, 1/2] . Several examples can be found in |fl8l . In this paper we refer to two of them repeat- 
edly. Those are the piecewise uniform Shishkin-mesh with <j)(t) = ItlnN and \j/(t) = N~ 2t , and 
the Bakhvalov-Shishkin mesh with 0(f) = -ln(l -2/(1 -AT 1 )) and y/(t) = 1 -2f(l -A^ 1 ). 
Besides above properties, we assume the function also to fulfil 



max (b'(t)<CN and min d 

fe[0,l/2] " i=L...,N/2 



N 



</> 



1 



N 



>CN 



-1 



We denote by Xij = [x{-i,Xi] x \yj-i,yj] a specific element and by x a generic mesh rectangle. 
Note that the mesh cells are assumed to be closed. Let /z,- :—xi — x/_i, kj :—yj —yj-i be the 
dimensions of Ty and 

h := max /?,-, k : = max k; and h m j„:= min h{. 
i=l,...,N/2 j=l,-,N/4 i=l,...,N/2 

Note that it holds, HH eq. (3.5)] 

h i <CeN~ 1 man\}ii , \e Px ^ ae \ i=l,...,N/2,xe[xt- h Xi], (2.3) 



max \w(t) 

te [0,1/2] 



and similarly for kj, where here and further-on max | \j/\ 

Let for a fixed polynomial degree p > 2 the finite element space be given by 

V N = {v e flj(Q) : v| T e Q p {%) Vt e T N } . 

Assumption 2.1. The solution u of (11.11) can be decomposed as 

U = V + W\ + W2 + Wl2, 

where for all xje [0, 1] and < i + j < p + 1 the pointwise estimates 

d i+ J Wl 



d'+Jv 

d'+J 
dx'dy 



<C, 



i*,y) 



<Ce- i e-M* 



dx'dy J 
and for i + j = p + 2 f/ze L 2 -norm bounds 

l|aja/v|| 0) n<c, 

||<?i<?>2||o,n<Ce^/ 2 + 1 /4 



dx { dyj 
^fcy) < Ce-ll 2 ( e ->'/v^ + e -(i->0/v^ 



d l+j wn , . 
r(x,y) 



(2.4) 



l^/wiHo^Ce-m/ 2 , 
|^/wi 2 || 0) fl<Ce-^/2)+3/4 



(2.5) 



/zo/J. Here w\ covers the exponential boundary layer, w 2 the characteristic boundary layers, 
w>i2 the corner layers, and v is the regular part. 

Remark 2.2. In / U21I731/ Kellogg and Stynes proved the validity of Assumption \2. 1 [ for constant 
functions b, c under certain compatibility and smoothness conditions on f. 

3 Interpolation 

We define two different interpolation operators. The first one tHEESl is the vertex-edge-cell 
interpolation operator ft p : C{x) — > Q p (x), defined locally on the reference element t : = [— 1, l] 2 
by 



{ft p v-v){di)=0, i'=l,..., 4, 



(3.1a) 
(3.1b) 



(7t p v-v)q = 0, i=l,...,4, qeP p - 2 (ei)i 

(7t p v-v)q = 0, qeQ p - 2 (t), (3.1c) 



where dj are the vertices and e,- the edges of t. Using the bijective reference mapping F T : f — > X, 
this operator can be extended to the global interpolation operator n p : C(Q) — > V N by 



-l 



-N 



(<v)| T := (ft p (voF z )) oF-\ Vt eT\ve C(Q) 



(3.1d) 



The second interpolation operator is of Lagrange-type. Let — 1 = to < t\ < ■ ■ ■ < t p — 1, be the 
zeros of 

(\-t 2 )L p (t)=0, re [-1,1], (3.2a) 

where L p is the Legendre polynomial of degree p. These points are also used in the GauB- 
Lobatto quadrature rule of approximation order 2p — 1. Therefore, we refer to them as GauB- 
Lobatto points. In literature they are also named Jacobi points [fT4l as they are also the zeros of 
the orthogonal Jacobi-polynomials P p of order p. 
The operator I p : C(t) — > Q p (f) is then defined on the reference element f by point evaluations 

(I p v) (t h tj) = v(t h tj), i, j = 0, . . . , p. (3.2b) 

With an extension like above we obtain the global interpolation operator Ip : C(Q.) — > V . The 
interpolation error for both operators can be bounded according to EKSl using Assumption ^. II 

Theorem 3.1. Let o > p + 1. Then it holds for the solution u of (11.11) 

1 1 \u -Ip u\ 1 1 £ < C(h +K + N- 1 max \\/\) p and \\\u - rf u\\\ £ < C(h+lc + N^ max |l/|) p . 

Remark 3.2. Note that 

,, I 21nA^, Shishkin mesh, 
max|v/|<<^ 

I 2, Bakhvalov-Shishkin mesh, 

which shows the improvement of the bounds using graded meshes near the boundaries. 

Lemma 3.3. Let ft p and I p be the local interpolation operators into Q p (f) on the reference 
element x defined in (13.11) and (13.21) . respectively. Furthermore, let ft p +\ be defined similarly to 
(13.11) interpolating into the local space Q p +\(t). Then it holds 

ft p = I p fc p+ \. (3.3) 

Moreover, ifl p+ i is a Lagrange-interpolant into Q p +i(i) using the interpolation nodes (t*,t*j) 
where {t*},i — 0,...,p + l consists of the p + 1 Gaufi-Lobatto nodes {f,}, i = 0,...,p from 
(I3.2al) plus one arbitrary node t* +l G ( — 1,1), then it follows 

h = ft p i; +l . (3.4) 

Proof We extend an idea given in ifTOll . For any function v EC(i) holds 

ftpV e Q p (x) and I p v e Q p (t). 

To prove the equivalence (13.31 ) we only have to show that I p ft p+ i shares the same degrees of 
freedom as 7t p . The definitions (13.11) and (13.21) imply 

(fc p v)(di) = (I p 7t p+ iv)(di) = v(di), i = 1, . . . ,4, (3.5a) 



where a,- are the vertices of t. Furthermore, it holds for q e P p _2(e,-) 

{v-I p 7t p+ iv)q= (v-l p 7t p+ iv-(v-7t p+ iv))q= (ftp+iV - I p (7C p+ iv))q -. V, 

where V E Pi p -\{ei) and e, is any of the four edges of t. Now, the integral can be rewritten 
using a quadrature rule that is exact for polynomials of order 2p — 1. We use the GauB-Lobatto 
rule and obtain 

f p 

/ (v-IpfCp + iv)q=Y,Wj((7tp + iv-Ip(7t p+ iv))q)\i i (tj)=0 : z'=l,...,4, (3.5b) 

Jit y= 

where {wj} are the weights of the quadrature rule. The last equality comes from (|3.2bl) . In a 
similar fashion it follows for q e Qp-iiX) 

(v - I p 7t p+l v)q = £ Wij((ftp + iv-Ip(fCp +l v))q)(t il tj)=0 1 (3.5c) 



due to the integrand being a polynomial in Q2 P -i (t), the GauB-Lobatto rule on the rectangle t 

and again (I3.2bl) . Comparing (13.51) to (13.11 ) one concludes (13.31) . 

The second equivalence (|3.4I) can be concluded easily by the first one: 

K plp+i = h K p+\Ip+i = Vp+i = h> 

where we use the property of I* +1 and ft p+ i being projections into Q p +i(t) in the second 
step. The last equality holds because I p uses a subset of interpolation nodes of I* +l in its 
definition. □ 



4 Supercloseness Analysis 

Let us now come to the numerical method. We define the Galerkin bilinear form by 

acai{v,w) :=£(Vv,Vw) + (cv-£v x ,w) 
and a stabilisation bilinear- form of the streamline-diffusion method IfTTTl by 

a s tab(v,w) := J^ 5 T (eAv + bv x -cv 1 bw x ) Tl 

T(ET N 

where the parameters <5 T > are user chosen and influence both stability and convergence. They 
are taken constant in each sub-domain of Q., i.e. d r = 5, ; - for any T C Cly. 
The streamline-diffusion bilinear-form is then defined as 

«5o(v, w) := a Ga /(v, w) + a stab (v, w) 
and the streamline-diffusion formulation of (11.11 ) is given by 



Find u N G V N such that 

as D {u\v N )=f SD {v N ), Mv N eV N (4.1) 

where 

/ffl>(v):=(/,v)- £ MfMx)v 

ter N 

Theorem 4.1. For Q p -elements, o > p + 1 and under the restrictions on the stabilisation pa- 
rameters 

8 n =CN~\ b\i <Cmax{l,e" 1 / 2 (A^" 1 max|v/|)}(iV" 1 max|v/ / |) 2 , 5 12 = b\ 2 = 

/zo/ds ?/ze estimate 

|||^a- M A '||| e <C(fi + ^ + ^V" 1 max|i/|) p+1 / 2 ln^V. 

Proof. In (3]| this result is given for the standard Shishkin-mesh. Together with techniques for 
S-type meshes, see e.g. |7l[T8||, the desired bound follows. Note that the additional logarithmic 
factor is caused by the estimation of the convective term inside the characteristic layers. □ 

Remark 4.2. By a close inspection of the proof of Theorem \4J\ in [3, Theorem 13] the sharper 
result 



IK 



^K-i^HI <C{h+K+N- 1 max\\i/\) p+l / 2 



can be given under the modified restriction 

52i <C(N^max\\/\) 2 . 
Corollary 4.3. Combining Theorems \3. 1 1 and \4. l\ yields the convergence result 

\\\u - u N \\\ £ < C(h+lc + N^ max\\/\) p . 
To analyse the supercloseness behaviour of the GauB-Lobatto interpolation operator, consider 

\\\tfu-U N \\\ £ <\\\rfu-U N \\\ £ + \\ \lp{Xp+lU - U) - {Tt^ +l U - U)\ I | £ + I \\Ttp +l U - U \ | | £ . 

(4.2) 



Its first term can be estimated by the supercloseness result of Theorem 14. II and its last term by 
the interpolation error result of Theorem 13.11 adapted to the case of elements of order p + 1 . 
Thus, we only have to estimate the energy norm of 

R N u:=I N p {7t N p ^u-u)-{K N p+l u-u). 

We start with some basic estimates for R N u. 



Lemma 4.4. For any w E C{Xij) holds the stability estimate 

II^HIm^)^ c II w IIm^)- < 43a) 

For any w E H P+2 (Q.), 1 < t < p and Tjj E T N we have the anisotropic error estimates 



r=0 
t+\ 

(^w) x || 0j% <c£ iw +1 - r ^i +w #Hlo*,, 

r=0 



(4.3b) 

(4.3c) 



and analogously for || (i^w^Ho,^- 

Proof. The stability estimate (I4.3al) is a direct consequence of the stability of the interpolation 
operators /r and ^r, j in L TO . Their stability holds because all degrees of freedom are point- 
evaluations or integrals. 
For (I4.3bl ) we use anisotropic error estimates [HJ0H7] to obtain 



dN 



R' y M\Q,rn<C 



<c 



r N 



JV 



hi(^p+iW-w) x \\o,x iJ + \\kj(7t 1 p+1 w-w)y\\o^ j 

rf+l t+l 

II",- k j d x ^^110,^ + 2^11"! k j d x d y W h,-Cij 

7=0 r=0 



which gives (I4.3bl) . For (I4.3cl ) we need additionally anisotropic estimates for the second order 
derivatives, see again 0] 



?N, 



(R n w) x \\o,r u <C 
<C 



t t 

||/j. fyd, OyWHo^+^ll"; */ »* «y w llo,^ 



(4.4) 



Lr=0 r=0 

which gives (I4.3cl) . D 

Theorem 4.5. Let us assume k <CN~ 1 / 4 , o > p+l and N~ 1 (max\\ir'\) 2 <C. For R N defined 
above holds 

\\\R N u\\\ £ <C(N-( a - l / 2 \l+h mesh ) + (h + k- + N- 1 m a x\x/\y +l ) 

where h mesh := eN~ l (\nN) 1 / 2 /h min . 

Remark 4.6. For the mesh specific value hmesh holds on a Shishkin mesh 

d(lniV)- 1 / 2 < h mesh < C 2 (\nN)-'l 2 , 
on a Bakhvalov-Shishkin mesh 

Ci(lniV) 1 / 2 <fc Me5/t <C 2 (lnA0 1/2 , 
and on a general S-type mesh 

d(ln^V)- 1 / 2 < h mesh < C 2 (\nN) 1 / 2 . 



Proof of Theorem 14.51 Let us start with the L2-norm estimate. We use the solution decomposi- 
tion of Assumption ^. 1 l and start with the regular part v. By (I4.3bl) with t = p— 1 > 1 we obtain 

ll**v|| ,o < C{h+k+N- l y+ l . (4.5a) 

Estimate (|4.3bl) can also be used to bound w\ in Q12 U Q22, where h\ can be estimated by (12.31) . 

^C^ + ArVaxlv/l) 2 ^ 1 ), (4.5b) 

where we used (p + l)/<7 — 1 < 0. In Qn U ^21 we use the stability (I4.3al) to obtain 

II^wiIIoaiu^! <C||wi|| M n llUn2l) <CAT ff . (4.5c) 

The other two layer terms can be estimated similarly and combining these results proves the 
L2-estimates. For the /^-component we use (I4.3cl) with t = p and its counterpart for the y- 
derivative for the regular solution component v to obtain 

||V2? iV v||o ) fl<C(/j+£+Ar 1 y H - 1 . (4.6a) 

Similarly, using the type of analysis as above, we show 

||V^w 1 || 0A2 ui2 22 <Ce- 1 / 2 (^ + iV- 1 max|v/|)P +1 , (4.6b) 

||Vi^W2||o,n M ufl 22 <C6- 1 / 4 (ft+Ar 1 max| V / |) p+1 , (4.6c) 

\\VR N w l2 \\o,£i 22 <C£- 1/4 (N- l m a x\\ir > \)P +l . (4.6d) 

On the other domains we use the decay of the layer terms. We show the analysis exemplary for 
the three terms yielding the largest bounds and invoking the most assumptions. Let us start with 
(R N wu)x in O-n. By (I4.3cl) with t — 1 we obtain 

< or^Ar^ArVaxli/DV 3 / 4 ]. 

On the other hand, a triangle and an inverse inequality give 

||(^VVi2)x||0,a 12 <C[^(||/^(^ +1 VVi2-VVi2)||0A2 + ll^+l W 12|l0A2) + ll( M; 12)x||0A2]- 



While the last term of the right-hand side can be estimated directly, we use an idea from [|20l 
incorporating the stability of /j; and ^ +1 for the other two terms. 

Pp (<+1^12 - W12) II0A2 + II ^+l w 12 II0A2 
< C ff (/X+1^12) 2 + ( 7 >12) 2 + «+lW 12 ) 2 
N „,- 3JV/4 v . 

<C £ / E / ( w 1 2(^-1 >3 ; ;-i)+ w 12(^-1 ^j))dy dx 

i=N /2+1 Jx '- 1 j=N /4+1 "'XM 

j=JV/4+l ^j-i V / 



<c 



N 






<c 



e £ dx+ e £ dx 
Jo 



j/v/4+1 _i^£ 



e 8 



1/2 



dy + 



1/2 _ 2v 



£ £ 



I 2 



<C[h + e] 



N~ l + e 



1/2 



2^/4 J %y 

N- 2a <Ce(e 1/2 +N- l )N- 2 °. 



dy 



Here we have used the symmetry of the pointwise bound of w\2 w.r.t. y. Thus, a second bound 
for [R W12)* in Q12 holds: 



fJ^w 



12jx||0,f2 I2 



<c 



h mesh £- l N(\nN)- l/1 N- a (£ l/ * +N- l/1 )£ 111 + £- l/4 N-° 
<C£- l l 2 N- G [£ l l* + h mesh {£ l l*N + N l l 2 )}. 

Combining these two estimates we obtain 

||(^w 1 2)J|oA2<Ce- 1/2 ^- a min{(^- 1 max|v/|) 2 e - 3 / 4 ,/ ?mei , /? (e 1 /4 A / +A /i/2 )} 

< Ce" 1 / 2 iV-( cr - 1 / 2 )(l +h mesh ), (4.6e) 



where we used iV _1 (max|y/|) 2 < C in estimating the minimum. The second term we want to 
look at is (R N W2)x in ^12- This one highlights in cancelling the logarithmic term, why h mes h is 
defined as it is. We obtain the two estimates 

||( J RV)x||oA2<Ce" 1/2 ^- f7 [(^+^- 1 ) 2 e- 1 / 4 ], 
||(tfV)*||oA 2 < Ct/w^MlnA^ 

<Ce-^ 2 N- a h mesh [e 1 ^N+N^ 2 ] 
and therefore 

||( J RV)x||oA2<Ce- 1/2 ^- a min{(^ + ^- 1 ) 2 e- 1 / 4 ,/ Jm , s ,[e 1 /^ + ^ 1 /2]} 

< C£- l / 2 N-^- l / 2 \l+h mesh ). (4.6f) 

As a third term we look at (R N wi) x in ^21 ■ A similar analysis as above gives for£ < CN~ l ' A 

\\(R N w l ) x \\ , sl2l <C£- l / 2 N- a mm{£- 2 ]c 2 1 l + (£N) 1 / 2 + £N} 

< C£-^ 2 N-°(\ + (M) 2 / 3 ) < C£-^ 2 N-^- l / 2 \ (4.6g) 



All other terms remaining can be estimated by similar steps. Combining the results gives finally 
the statement of the theorem. □ 

Remark 4.7. Following the proof of Theorem 14.51 and neglecting the restriction on It, one can 



still show with ( 4.6g ) for o > p + 1 

\\\R N u\\\ <C(iV-( fT - 2 / 3) + (ft+£+Ar 1 max|y/|)P +1 ). 



Theorem 4.8. Let o > p + 2. Then it holds for the streamline-diffusion solution u under the 
restrictions on the stabilisation parameters given in Theorem \4.1\ 

|||/pW-^||| £ <C(fi + £ + ArVax|i/|) p+1 / 2 lnAr. 

Proof. Consider again (14.21 ) 

lll/^« — W^lll < III^M-M^HI + |||/^(^, ,U-u) — {l£< iW — w)||| + Hl^i iW — will . 
Ill p 1 1 1 e — IMP 1 1 1 e 1 1 1 P ^ P+ l ' K P+ l ' 1 1 1 e 1 1 1 P+ l 1 1 1 e 

Theorem 14. II gives under conditions on the stabilisation parameters and o > p + 1 



\*p- 



- N u-u N \\\ £ <C(h + lc + N- l max\\j/\)P +l / 2 \nN, 
Remark 14771 gives for o > p + 5/3 

1 1 \lp(7$ +l u -u)- (rf +i u -u) 1 1 1 < C(h + K + N- 1 max \tf\) p+ \ 
and Theorem |3T| yields for o > p + 2 

|||7^ +lM - M ||| £ <C(/J+£ + Ar 1 max|v/ / |) p+1 . 

Combining the three estimates completes the proof. □ 

5 Superconvergence by Postprocessing 



By utilising the supercloseness results of Theorems 14. 1 I and l4.8l we can construct postprocessing 
operators. They improve our numerical solution with little additional computational effort to 
higher convergence order. 

Suppose Af is divisible by 8. We construct a coarser macro mesh f N ' 2 composed of macro rect- 
angles M, each consisting of four rectangles of T N . The construction of these macro elements 
M is done such that the union on them covers Q. and none of them crosses the transition lines 
at x = X x and at y = Xy or y = 1 — Ay, see Figure |2l Remark that in general f N ' 2 ^ j N l 2 due to 
different transition points X x and X y and the mesh generating function . 
We now define postprocessing operators locally for one macro element M G T N ' 2 . The first one 
was presented in Id in fBTTl and is a modification of an operator given in lfT5l . 
Let v be the linearly mapped function v from any interval [jc/_i,x ( ' + i] onto the reference interval 
[—1,1]. Note that xi is not necessarily mapped onto 0, but to a value aG (-1,1). In [2] 
the following condition on the underlying mesh is given, that guarantees the non-degenerate 



\-Xy 



Ax 



Figure 2: Macroelements M of f N < 2 constructed from T 



N 



behaviour of the macro elements and of the operators defined on it: There exists a constant 
q > 1 independent of N and e such that 



max{hi,h i+i } 
mm{hi,h i+l } 
max{kj,kj + i} 
mm{kj,k j+ i} 



< 



for all i= I,..., N/2-l, 
<q, forall j= l,...,N/4-l and j = 3N/4+l,.. .,N -1. 



A Shishkin mesh has q = 1 while a Bakhvalov-Shishkin mesh has g = ln(3)/ln(5/3). For many 
more S-type meshes this condition holds and we assume it further-on. 
Define the reference operator P vec : C[— 1 , 1] — > P p +\ [—1,1] by 



for p = 2: 
while for p > 3: 



-l 

a 



P vec v(-l) = V(x/_l), PvecV{a) = v(xi), PvecV(l) = v(x i+ l), 

(P vec v-v)=0, 



(P ve ct-f) = 0, / (P vec v-v) = 



(P vec v-v)p = 0, p6P«_ 2 [-l,l]\: 



By using the reference mapping and the tensor product structure we obtain the full postprocess- 
ing operator P vec .M '■ C(M) — > Q p +\ (M) on each macro element. Then, this piecewise projection 
is extended to a global, continuous function by setting 

(PP+'v) (x,y) := (P vec , M v) (x,y) for (x,y) e M. 

The second postprocessing operator is defined by using only point evaluations. Let {(xi,yj)}, 
i,7 = 0,..., 2p denote the ordered sample of GauB-Lobatto points of the four rectangles that M 
consists of. Let Pql.m '■ C(M) — > Q p +i (M) denote the projection/interpolation operator fulfilling 



PcL,M v ( x i^j) = v (*iJj), hJ = Q, 1,3, 5,..., 2/7-3,2p-l,2p. 



Then, this piecewise projection is extended to a global, continuous function by setting 

(pP+'v) (x,y) := {Pglmv) (x,y) for (x,y) G M. 
Lemma 5.1. For the postprocessing operators Pq L and P$ ec defined above holds 



P P G l l I N p v = P p G + L \ 



P P ot l v N 



<C v 



N\ 



P P+l K N V = P P+1 V 
r vec ll p v r vec v i 

\p p + x v N \\\ <cllliAIII 

\ r vec v lie— Mil lie' 



forallveC(Q.), 

for all v N eV N . 



Let u be the solution of (11.11) . Assumption \2. 1 1 be true and o > p + 2. Then it holds 



GL 



+ 1 1 \PP+ I u -u\\\< C(h + K + N- 



max 



\V\) 



p+1 



(5.1a) 
(5.1b) 

(5.1c) 

Proof. The consistency (I5.1al) is a direct consequence of the definitions of Pq L and P p ec ■ 
The stability (I5.1bl) can be shown for both operators similarly. Therefore, let Pm be any of the 
local operators Pql,m and P vec .M- The stability (I5.1bl) then follows by showing ||-Pmv^||o.a/ — 
C||\A|| ^ and \Pmv*\i,m < C|v%,m for any M G f N/2 and any v^ G V N (M). The operator P M 
is a linear operator from the finite dimensional space V N (M) into the finite dimensional space 
Q p+ i(M). Thus it is continuous and with v N H> ||v^||o,m being a norm in both spaces we obtain 



\\Pmv 



.N\ 



|0,M 



N r- \/N I 



<C\\v n \\qm, for all v^GV^M). 



Similarly v H- |v |im is a norm on the quotient space V (M) \Mand v i-» ||v ||i,m is a norm 
on Q p+ i(M). Therefore, 



\Pmv n \i,m < \\Pmv n \\i,m < CK|i )M for all v N E V N (M)\R. 

Finally, the interpolation error ()5 . let follows by Assumption 12. II and standard anisotropic esti- 
mates for interpolation [THE). d 

Theorem 5.2. Let o > p + 2. Then it holds for the streamline-diffusion solution u N under the 
restrictions on the stabilisation parameters given in Theorem \4.1\ 

+ 1 1 \u -P p + l u N \ 1 1 £ < C{h+R + N- 1 max | v/|) p+1 / 2 ln^V. 
Proof. Using the consistency and stability of Pq L we obtain 



u-P^u" 



u-p£W 



< 



< 



U- P GL U 

u — Pq L u 



+ 



Pof^-Pot^ 



i /-> 1 1 1 tN N\\\ 

+ C \\\r„u — u \\\ , 
III p I Me 



Similarly one can show 



_pP+l,.N\ 



< \\\u-PP+ l u\\\ p + C\\\^u-u 



N\ 



Now the statement follows by (I5.1cl) and the supercloseness results of Theorems 14.11 and 14.81 

□ 



Table 1: Convergence and closeness errors for ^-elements on a Shishkin mesh for e = 10 

iV 



-6 



with corresponding rates r^ 



N 


W — U N 


1 1 \F. 


fl^u — u 


AMI 

\\\F. 


\\\l$u-u' 


VIM 
\\\f 


J^u — u 


N \ 

\\\f 


8 


2.270e-02 


2.60 


8.259e-03 


2.55 


9.587e-03 


2.72 


9.811e-03 


2.68 


16 


7.926e-03 


2.79 


2.940e-03 


2.71 


3.184e-03 


2.81 


3.309e-03 


2.79 


32 


2.141e-03 


2.96 


8.233e-04 


3.28 


8.515e-04 


3.31 


8.923e-04 


3.24 


64 


4.723e-04 


3.04 


1.544e-04 


4.25 


1.572e-04 


4.24 


1.707e-04 


3.83 


128 


9.156e-05 


3.01 


1.563e-05 


4.53 


1.600e-05 


4.50 


2.169e-05 


3.46 


256 


1.699e-05 




1.240e-06 




1.292e-06 




3.122e-06 





Table 2: Convergence and closeness errors for (^-elements on a Bakhvalov-Shishkin mesh for 
£ = 10~ 6 with corresponding rates r^ 



N 


w — u N 


1 1 \f 


K^U — U 


AM 
\\\f 


\l^u — u' 


v 
1 1 If 


\j^u — u 


N \\ 

\\\f 


8 


4.104e-03 


2.74 


1.673e-03 


3.07 


1.774e-03 


3.11 


1.828e-03 


3.05 


16 


6.145e-04 


2.92 


1.986e-04 


4.22 


2.055e-04 


4.18 


2.207e-04 


3.67 


32 


8.106e-05 


2.94 


1.065e-05 


4.33 


1.137e-05 


4.25 


1.729e-05 


3.21 


64 


1.055e-05 


2.97 


5.309e-07 


4.17 


5.960e-07 


4.11 


1.868e-06 


3.02 


128 


1.349e-06 


2.98 


2.950e-08 


4.08 


3.440e-08 


4.05 


2.310e-07 


2.99 


256 


1.707e-07 




1.742e-09 




2.077e-09 




2.902e-08 





2 



6 Numerical Example 

Let us consider the singularly perturbed problem given by 

-£Au-(2-x)u x + 3/2u = f, inH = (0,1 

u = 0, on dQ., 

with a constructed right-hand side, such that 

e -x/e_ e -l/e\ ( l _ e -y/e^ ) ( l _ e -(l-y)/eVi) 



u(x,y) — cos(x^:/2) — ■ 



1-e-Ve / l-e-!/£ 1/2 



is the exact solution. 

The following calculations were done in Matlab and the linear systems solved by its direct "\"- 

solver. We fix the polynomial degree to p = 3 and the parameter for the Shishkin mesh to a = 

p + 2 — 5. Moreover, we set £ = 10 , sufficiently small to generate the sharp boundary layers 

we are interested in. Note that additional computations were done with different polynomial 

degrees p and varied perturbation parameters £ supporting the same conclusions. 

In the following tables, the experimental rates of of convergence for given measured errors e^ 

are calculated by 

s _ ln(e N /e 2N ) B _ ln(e N /e 2N ) 

PN ~\n(2\n(N)/\n(2N)) , Pn ~ ln(2) ' 



Table 3: Superconvergence errors for Qt, -elements on a Shishkin mesh for e = 10 6 with cor- 
responding rates rfj 



N 


1 1 U P N 1 

1 1 u r vec l 


'111, 


1 1 a P N i 

\\\U — l/^I " 
1 UL 


'111, 


8 


2.940e-02 


2.79 


2.883e-02 


2.77 


16 


9.471e-03 


3.16 


9.362e-03 


3.21 


32 


2.142e-03 


3.54 


2.069e-03 


3.59 


64 


3.511e-04 


4.03 


3.298e-04 


3.96 


128 


4.001e-05 


4.11 


3.892e-05 


4.06 


256 


4.005e-06 




4.005e-06 





Table 4: Superconvergence errors for (^-elements on a Bakhvalov-Shishkin mesh for e = 10" 
with corresponding rates rfj 

N IIIm-P^ u n \\\ \\\u — P n ij n \\\ 

I I I v L L ! I I r, I I I \J * ' I I I fc 

16 8.965e-04 3.99 9.771e-04 3.71 

32 5.633e-05 4.07 6.141e-05 3.99 

8 1.071e-02 3.58 1.282e-02 4.10 

64 3.343e-06 4.06 3.587e-06 4.09 

128 2.003e-07 4.04 2.110e-07 4.06 

256 1.217e-08 1.266e-08 



assuming epj = C(N~ l lnN) p N on a Shishkin mesh and e^ = CN~ p N on a Bakhvalov-Shishkin 
mesh. 

The numerical method is the SDFEM given in (14.11) with stabilisation parameters according to 
the upper bounds in Theorem 14. II where C is set to 1. 

Tables [Hand [2] present the convergence and closeness results. We observe third order conver- 
gence of the numerical method as predicted by Corollary 14.31 Moreover, the results on a graded 
mesh like the Bakhvalov-Shishkin mesh are much better compared with the piecewise equidis- 
tant Shishkin mesh. There are two orders of magnitude difference in the final line for the energy 
error. 

The closeness results for the vertex-edge-cell interpolant IK and the GauB-Lobatto interpolant 
1^ as well as a pointwise interpolation operator J^ using equidistantly spaced interpolation 
points are also given. Supercloseness of a better order than p + 1/2 = 3.5 can clearly be seen 
for 7E3 and 1^, whereas for J^ the rate is not as high. On the Bakhvalov-Shishkin mesh, we 
observe a clear order 4 for the first two operators and an order 3 for the equidistant interpolant. 
Tables [3] and [4] now show the results of the postprocessed numerical solutions. We observe 
for both postprocessing operators on both meshes convergence rates of order p + 1 =4, which 
compared with the convergence results presented in Tables [Hand [2] is an increase of a full order. 
Theorem 15.21 only predicted half an order increase. Thus, it seems that the supercloseness result 
of Theorem 14. II is not sharp. By improving this estimate, the improvement of Theorem 15.21 



follows immediately. 

By comparing the results of the two operators, we observe very little difference for larger values 
of N. Thus the new postprocessing operator using only point values has comparable conver- 
gence properties to the already existing one that uses integral values. 

Acknowledgement: The author would like to thank Lars Ludwig for carefully reading the 
manuscript and suggesting improvements to the readability of the paper. 
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